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ABSTRACT 

We present observations of ^^CO(l-O) in 17 Combined Array for Research in Millimeter 
Astronomy (CARMA) ATLAS^^ early-type galaxies (ETGs), obtained simultaneously with 
^^CO(l-O) observations. The in six ETGs is sufficiently bright to create images. In 

these 6 sources, we do not detect any significant radial gradient in the ^^CO/^^CO ratio be¬ 
tween the nucleus and the outlying molecular gas. Using the channel maps as 3D masks 
to stack the ^^CO emission, we are able to detect 15/17 galaxies to > Scr (and 12/17 to at least 
5 < t ) significance in a spatially integrated manner. Overall, ETGs show a wide distribution of 
^^CO/^^CO ratios, but Virgo cluster and group galaxies preferentially show a ^^CO/^^CO 
ratio about 2 times larger than field galaxies, although this could also be due to a mass depen¬ 
dence, or the CO spatial extent (i?co/A)- ETGs whose gas has a morphologically-settled 
appearance also show boosted ^^CO/^^CO ratios. We hypothesize that this variation could be 
caused by (i) the extra enrichment of gas from molecular reprocessing occurring in low-mass 
stars (boosting the abundance of ^^C to ^^C in the absence of external gas accretion), (ii) 
much higher pressure being exerted on the midplane gas (by the intracluster medium) in the 
cluster environment than in isolated galaxies, or (iii) all but the densest molecular gas clumps 
being stripped as the galaxies fall into the cluster. Eurther observations of ^^CO in dense en¬ 
vironments, particularly of spirals, as well as studies of other isotopologues, should be able to 
distinguish between these hypotheses. 

Key words: galaxies: clusters: general - galaxies: elliptical and lenticular, cD - galaxies: 
ISM 


1 INTRODUCTION 

While early-type galaxies (ETGs) wer e originally assume d to be 
mainly devoid of molecular gas (e.g.. ISower et aljll992h . it has 
become increasingly apparent that these “red and dead” early- 
type galaxies often contain non-negligible amounts of cold gas 

jWiklind & Rvdbecldl 19861:IPhillips et ^Jl 19871:Iwiklind fc Hentell 

* E-mail: kalatalo@caltech.edu 


I^: ISage & Wrob^ 1 19891: IWelch & Sag3 l2003l : ISage et all 
2007h . Large-scale searches for molecular gas in ETGs have indeed 
shown that a significant proport ion of these galaxies do contain 
a mo lecular gas reservoirs (e.g.. ICombes et ^l2007l : [Young et al.l 
I 2 OI ih . though many of the early searches targeted relative small 
unrepresentative samples of ETGs, rather than unbiased, volume- 
limited samples. 

The ATLAS®° survey jCapnellari et alj [ 2 OI lah provides a 
complete, unbiased sample of morphologically-selected ETGs in 
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Figure 1. (Left): Distribution of luminosities (gray ) among both ^®CO-detected ETGs (green) and this work (blue). The comparison sample is 

from the flux-limited survey of lYoung et al.l i201 ill . The sample represents a large percentage of ETGs detected at Lqq > 10^ ® Lq, though 

misses some of the faintest ETGs detected in ^^CO. (Right): Distribution of all galaxies with measured including ETGs (gray). The ^®CO-detected 
ETGs are distinguished (blue) along with galaxies part of this work (green). In this c ase, while the ETG s fall in the low Lqq end of the distribution, it appears 
that this is due to the fact that ETGs in general have lo wer Lgo (fur ther detailed in lYoung et al.ll20il1) . rather than the ^®CO-detected ETGs being biased. In 
both cases, the smallest value of Lco NGC 404 from lSag3 jl990ll . atLco~ 6 L 0 . 


a large local vol ume. All ATLAS^^ galaxies were searched for 
molecular gas by lYoung et al.l ( 1201 ih . with a measured detection 
rate of 22 ± 3%. The detection rate of molecular gas did not appear 
to be dependent on galaxy environment, with a similar fraction of 
group and cluster ETGs de tected as field ETGs. Follow-up efforts 
included CO imagi n g (e.g., Youn g et '^120081 : ICrocker et alJlOlll : 
lAlatalo et alj|2013l) . Davis et al. 1 201 ih found that despite the fact 
that Virgo galaxies are detected at a similar rate to field ETGs, the 
origin of the gas within Virgo galaxies appears to be exclusively in¬ 
ternal, while a large fraction of field ETGs acquire their molecular 
gas via accretion of gas-rich companions. 

*^^ 00 ( 1 - 0 ) has long been used as an optically thin tracer 
of ^^CO, and for the most part has been shown to be a faith¬ 
ful tr a cer of it l Young & Sandersiflosd : ISagelll99^ : ISage & Isbelll 


Il99ll: lAaltoetalJll995[) . The range of i®CO/^^CO ratios (here¬ 
after ^ 13 ) seen in past studies (mainly in spirals, starbursts and 
mergers) was 0.06-0.4, with ong oing major merg e rs populatin 
the lower end of the distribution dAalto et ^Il99lh . ISagel h99' 
studied the i®CO(l-0) properties of 3 ETGs (NGC404, NGC4710 
and NGC 5195) with the National Radio Astronomy Observatory 
(NRAO) 12m single-dish, and concluded that this small sample of 
SO galaxies seems to have a similar value of as spiral galax¬ 
ies on average. NGC 5195, known to have recently interacted with 
M51, has a much smal ler ,^13 than the oth e r ETGs, similar to w hat 
is seen in NGC 1266 dAlatalo et al.ll 201 H : ICrocker et_ay j20lj)^ 
pecul iar ETC with compact, dense molecular gas 1 Alatalo et al.l 
l2015h currently hosting an active galactic nucleus (AGN)-drive n 
massive molecular outflow dPavis et al.ll2012l : lAlatalo et alj|2014l) . 


Dense gas (i.e., gas that is traced by HCN and HCO'*' emis¬ 
sion) was obs erved in a much larger sample o f ETGs through 
the SAURON dSacon et al ]|200lh and ATLAS®° dCappellari et al.l 
l 2011 ah surveys, that targeted the strongest ^^CO-detected 
ATLAS^^ galaxies with the Institut de Radi oastronomie Mil- 
limetrique (IRAM) 30m single-dish telescope dKrips et al.l 120101 : 


ICrocker et alj |201 IPavis et ^l2013l) . In these works, there was 
some evidence that while ETGs show similar ratios of dense gas as 
spirals, the ratios may depend on environment, with Virgo galax¬ 
ies showing boost ed ,^13 and HCN/CO ratios compared to field 
galaxies. With the ICrocker et al.l d2012h sample, statistically sig¬ 
nificant correlations between ,^13 and absolute A'-band luminosity, 
stellar population age, molecular to atomic gas ratio, dust morphol¬ 
ogy and dust temperature were found. The correlation between ,^13 
and environment was not however statistically significant, perhaps 
due to the low numbers of galaxies in dense environments such as 
the Virgo cluster. 

Here we present the ^®CO(1-0) observations from the 
CARMA array. In Sj3] we discuss the simultaneous ^^CO(l-O) and 
^®CO(1-0) observations, and the techniques that were used to re¬ 
duce and analyse the ^®CO(1-0) data, including imaging strong 
detections and stacking all non-detections. In we describe our 
sample galaxies and their properties. In |j4l we discuss the implica¬ 
tions of the distributions and how our new detections inform 
upon what affects ,^13 in galaxies. In SjS] we summarize our con- 
cl usions. The ATLAS^^ sa mple properties are described in detail 
in ICappellari et al.l ( 1201 lalF , and details s pecifi c to the CARMA 
observations are reported in lAlatalo et al.l ( l2013h . We use the cos- 
mological paramete rs Hq = 70 km s“^, Dm = 0.3 and A = 0.7 
dSpergel et alj2007l) throughout. 


2 THE SAMPLE 

The CARMA sample has been extracted from the ATLAS®^ sam¬ 
ple of galaxies, which is a volume-limited sample (D<42Mpc) 
of massiv e (Mx <- 21 .5 mags) mor phologically-selected early-type 
galaxies dCappellari et alj boi lah . The subsample of galaxies 


^ http://www.purl.org/atlas3d 
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Figure 2. Integrated intensity map of the ^^CO(l-O) (blue contours) overlaid upon the momentO map (colourscale) for NGC3665, NGC4429, 

NGC4435, NGC4710, NGC4753 and NGC5866, the galaxies with strong (see ^3.U detections. The CARMA ^^CO(l-O) beam is shown at the 

bottom-left of each panel, indicating that the galaxies are well res olved, and a scale bar re presenting 10^' is shown on the bottom-rig ht. Galaxies known to be 
in the Virgo cluster include NGC 4429, NGC 4435 and NGC 4710 |Can£ellmT^t^ .1201 lah . NGC 4753 is well-known to be in a group jSteiman-Cameron et alj 
Il992h and NGC 3665 and NGC 5866 are within a galaxy association iGarcialll993ll . The^^CO levels are [0.1, 0.3, 0.5, 0.7, 0.9] times the maximum. 


Table 1. ^®CO(1-0) observations for the CARMA ATLAS®° survey 


Name 

^maj ^ ^min 

Velocity Range 


0 

0 

-P’CARMA ^ 
Faota 

.^13 

pi'" 

'^1200 

/ 60//100 


(") 

(km s-l) 

(km s“^) 

(Jy km s ^) 


(Jy km s“^) 


IC676 

4.57 X 4.00* 

1360-1520 

20 

6.52 ± 1.78 

0.96 ±0.3 

0.096 ±0.03 

67.9±5.8 

0.62“ 

IC719 

4.09 X 3.74 

1670-2080 

40 

3.45 ± 0.64 

- 

0.14 ±0.03 

23.8±2.9 

0.33 

NGC 2764 

4.14 X 3.79 

2570 - 2840 

30 

9.05 ± 1.24 

1.4 ±0.2 

0.11 ±0.02 

79.7±4.2 

0.51“ 

NGC 3607 

5.82 X 5.19 

680- 1210 

40 

8.57 ±0.69 

1.1 ±0.1 

0.16 ±0.02 

53.2±3.2 

0.79“ 

NGC 3619 

4.65 X 4.11 

1280 - 1770 

80 

< 1.25 

- 

< 0.083 

15.1±1.3 

0.21 

NGC 3626 

4.14 X 3.84 

1320-1670 

50 

3.89 ± 0.83 

- 

0.097 ±0.02 

40.1 ±3.4 

O.SO*” 

NGC 3665 

4.42 X 4.27 

1750-2400 

50 

20.0 ±0.80 

1.1 ± 0.07 

0.22 ±0.01 

91.2±4.4 

0.30“ 

NGC 4324 

4.93 X 4.04 

1520 -1770 

40 

3.06 ±0.50 

- 

0.15 ±0.03 

20.0±2.1 

0.21 

NGC 4429 

5.12 X 3.90 

830- 1370 

30 

20.5 ± 0.66 

- 

0.29 ±0.02 

70.4±3.2 

0.30 

NGC 4435 

4.50 X 4.19 

560 - 980 

30 

5.97 ±0.40 

- 

0.18 ±0.02 

32.6±2.6 

0.43 

NGC 4694 

4.28 X 3.39 

1150-1210 

30 

3.83 ±0.63 

2.3 ± 0.7 

0.23 ± 0.04 

17.0±0.8 

0.41“ 

NGC 4710 

4.37 X 4.07 

940-1310 

30 

51.5 ± 2.47 

2.3 ± 0.1 

0.15 ±0.007 

350.6±6.0 

0.39“ 

NGC 4753 

5.81 X 4.32 

919-1470 

50 

12.7 ± 1.05 

- 

0.18 ±0.02 

70.2±4.9 

0.29 

NGC 5379 

5.65 X 3.90 

1700-1880 

30 

1.83 ±0.42 

- 

0.080 ±0.02 

22.9±2.3 

0.21 

NGC 5866 

3.75 X 3.27 

500 - 1010 

30 

34.1 ± 2.59 

2.2 ± 0.2 

0.14 ±0.01 

250.9±8.4 

0.27“ 

PGC 029321 

4.10 X 4.01 

2760-2910 

30 

< 1.46 

- 

< 0.074 

19.6±1.7 

1.3 

PGC058114 

4.72 X 3.93 

1370-1640 

30 

5.68 ± 1.13 

2.1 ± 0.5 

0.077 ±0.02 

73.4±3.1 

0.78“ 


t Channel width of the '^^CO cube. 

t Does not include the 20% absolute millimeter flux calibr ation uncertainty. Errors were calculated as discussed in am 
^ 30m integrated flux densities jCrocker et al.ll2Q12h . 

^ integrated flux densities from the CO cubes described in lAlatalo et alj i2013l) . but derived with the same velocity resolution (and channels) as the 

cubes to match the resolution for our st acking analys is. _ 

*The map of IC 676 is uniform ly weighted in lAlatalo et al.l i2013h . whereas the other and maps in this work and are naturally weighted 
“/eo/ZlOO from lCrocker et alj j20T^. Others values are calculated from IRAS fluxes. 

^IRAS data were unavailable, so the flux ratio was calculated based on the Herschel /70//16O ratio (T. Bitsakis, private communication). 
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observed with CARMA ha ve been sele c ted a mong those de¬ 
tected in the line by lYoung et alj bOl ih . who observed 

the complete ATLAS®° sample. CARMA then imaged all pre¬ 
viously un-imag ed CO-detected E TGs above a flux limit of 
18.5 Jy km s“^ jAlatalo et al.ll2013h . Figure [TJi shows the distri¬ 
butions of luminosities for ^^CO-detected ETGs from lSag^ 
(Il990h : ICrocker et al.l ( 1201 2h and this work com pared to the un¬ 
biased, flux-limited survey of I Young et al.l ( 1201 ih . indicating that 
^®CO-detected ETGs tend to populate the high Leo end of the dis¬ 
tribution. Given that CARMA obs erved the brightest ^^CO ETG 
detections from lYoung et al.l ( 1201 ih . the fact that ^®CO detections 
tend to populate the bright portion of the ^^CO distribution is un¬ 
surprising. 

Eigure[T] shows that when we draw compa risons to the avail¬ 
able ^ ^CO-detected galaxies in the literatu re dYoung & SandersI 
Il986l: ISage & Isbelllll99ll: lAalto et al.lll995h . that are either late- 
type galaxies (LTGs), starbursts or mergers, the ^®CO-detected 
ETGs fall to the lower end of the distribution of ^^CO luminosi¬ 
ties, consistent with wh at is found in ^^CO-detected ETGs in gen¬ 
eral dYoung et al.ll201 ih . So our ^^CO ETG sample has properties 
entirely consistent with a population of ^^CO-bright ETGs. 

It also appears clear that the galaxies added in this work do 
not fundamentally differ in their ^^C O properties from t he galax¬ 
ies considered in previous ETG work dYoung et alj|201 ih . and ap¬ 
pear to be a robust addition to the ^®CO-detected ETGs. If we 
look specifically at the representation of Virgo galaxies, we find 
that 7/26 (26±8%) are Virgo members within the ^®CO-detected 
sam ple, compared to 12/56 (21±5%) of the ^^CO-detected ETGs 
from lYoung et alj d201lh . Thus it does not appear that Virgo mem¬ 
bers are over- or under-represented in our sample. ^^CO-detected 
ETGs were also selected from this subsample, using identical se¬ 
lection criteria to Virgo and thus are also neither over- nor under¬ 
represented in the sample. 

Though the 7 Virgo members have comparable average prop¬ 
erties to the non-Virgo ^®CO-detected ETGs (within the scat¬ 
ter), including stellar mass (< MK.vugo >= —23.6 ± 0.9 vs. 
< ALk, non-Virgo >= —22.9 ± 1.1) and molecular gas extent 
(< Rco /Re. viTeo >= 0.45±0 .45 vs. < Rco/Re ,non —Virgo ^ — 
0.67 ± 0.43: IPavis et 'ZI I2013I) . they do skew more massive with 
more compact molecular discs. This is likely due to the small num¬ 
ber of Virgo galaxies within this sample, but the scatter is large 
due to small number statistics. Eurther ^^CO studies of more Virgo 
members should reduce the scatter. 


3 OBSERVATIONS AND RESULTS 

^^CO(l-O) observations were obtained for 17 of the 30 original 
CARMA ^^CO(l-O) galaxies. The galaxies missing ^®CO(1-0) 
data are those imaged before Spring 2009, when upgrades to the 
CARMA receivers allowed for simultaneous imaging of both lines. 
The total available bandwidth from CARMA in each line varies 
as well depending on the semester the source was observed, due 
to upgrades to the c orrelator. Tota l bandw idths for each object are 
listed in Table 2 of I Alatalo et al J ( l2013l) . All data reduction and 
calibratio n were done s i multa neously with the ^^CO imaging pub¬ 
lished in lAlatalo et al.l ( l2013tl . Table m presents the galaxies ob¬ 
served in ^®CO, along with ^®CO-specific observing parameters 
(including beam size, bandwidth and channel widths). The calibra- 
tors used, and total observing time for each source, can be found in 
lAlatalo etal]j2013l) . 


Table 2. Resolved ^13 values 


Name '^i3,nuc -^ 13 ,ext 


NGC 3665 

0.24 

± 

0.01 

0.29 

± 

0.01 

NGC 4429 

0.29 

± 

0.01 

0.41 

± 

0.02 

NGC 4435 

0.35 

± 

0.01 

0.34 

± 

0.01 

NGC 4710 

0.18 

± 

0.01 

0.17 

± 

0.01 

NGC 4753 

0.15 

± 

0.03 

0.10 

± 

0.03 

NGC 5866 

0.06 

± 

0.02 

0.08 

± 

0.02 


3.1 Strong detections: maps 

The ^®CO(1-0) emission was strong enough (emission was de¬ 
tected with signal-to-noise ratio >3) to map 6 of the 17 observed 
galaxies. M oment maps were c onstructed based on the methods 
described in lAlatalo et alj ( l2013h . with a full velocity width deter¬ 
mined by the ^^CO(l-O) line. Eigure|2]shows the six mapped galax¬ 
ies: NGC 3665, NGC4429, NGC4435, NGC4710, NGC4753 and 
NGC 5866. To test the agreement b etween the CARMA detections 
and the single dish detections from ICrocker et all (l2012h . we con¬ 
volved the images of the three overlapping galaxies with strong 
i^CO detections: NGC 3665, NGC 4710 and NGC 5866, to the size 
of the IRAM 30m beam, and compared ,^ 13 , finding them t o be i n 
agreement within errors with ,^13 reported in ICrocker et al.l ( l 2012 h . 

For these six galaxies, we attempted to investigate whether 
,^13 varies with position. We first regridded both the ^®CO and 
^^CO maps so that the area of one square pixel equalled the area 
of the ^®CO beam, this using the Interactive Data Language (IDL) 
Astronomer’s Library routine hcongricfl with bilinear interpo¬ 
lation of all points. We then measured both a nuclear (from the 
regridded pixel corresponding to the 2 Micron All Sky Survey Kg 
band-determined nuclear position: Iskrutskie et alj2006h and a non¬ 
nuclear .^13 value (from the remaining ^^CO-detected pixels). To 
evaluate the uncertainty on the .^13 measurements, we first evalu¬ 
ated the uncertainty on each integrated flux by measuring the root 
mean square noise (RMS) in the regridded ^^CO and '^^CO inte¬ 
grated intensity maps, respectively, taking the standard deviation of 
the pixels outside of the region where '^^CO is detected, and multi¬ 
plying by v/iV, where N is the number of pixels in each region (one 
for nuclear regions; the total number of detected off-nuclear pixels 
for the non-nuclear regions). These M 13 gradient measurements are 
listed in Table[2]and discussed in Section|4l 

3.2 Weak ^®CO detection: ^®CO spectra 

To extract ^^CO integrated flux densities for the fainter galaxies, 
where ^^CO is not detected in individual beams, we utilize our prior 
knowledge of the molecular gas distribution from the brighter ^^CO 
emission. This technique is also used to extract integrated flux den¬ 
sities for the strong '^^CO detections. We used a smoothed ^^CO 
datacube (created to correspond both to the weighting scheme as 
well as velocity properties of the ^®CO maps) for each galaxy 
as an priori dataset, and for each channel determined the two- 
dimensional region over which the line is detected with a 

signal-to-noise ratio abov e 1.5. This 3D regi on then represents a 
detection mask tsee lSchruba et al .I2OIII for details of a simi¬ 
lar method). For each galaxy, we then sum the datacube over 


^ http://idlastro.gsfc.nasa.gov/ftp/pro/astrom/hcongrid.pro 
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Figure 3. Stacked ^®CO(1-0) spectra of the 15/17 detected galaxies in the CARMA ATLAS®® ®®CO survey. The dark blue shaded regions represent ztlcr 
uncertainties, using the noise calculation laid out in j]L2J w ith the light blue representing the integrated signal within the 2D ^'^CO(l-O) mask of each channel, 
compared with the ^-^CO/l-O) spectra from lAlatalo et in red. The un-shaded spectra are those channels where the number of pixels in the 

mask fall to fewer than total number of pixels of the beam. In these channels, the spectrum is calculated by summing pixels in a contiguous region of the 
channel map the size of the beam area, and can be considered the baseline scatter for the stacked spectra, with con'esponding error bars in the top left of 

each panel. The scale on the left of each panel is the ^^CO(l-O) line flux, while the scale on the right of each panel is the ^^CO(i-O) line flux. The velocities 
are the optically-defined , local standard of rest ('Uopt,LSR) velocities of the molecul ar gas. Total bandw idths from the CARMA observations are listed in Table 
2 of lAlatalo et al.l i2Q13l) . In all cases except NGC 3665, the continuum (reported in lAlatalo et al.l2Q'l3h is insubstantial compared to the stacked ^^CO spectra 
(and has been removed from all line data). A dashed black line shows the continuum detection level of 8.98 mJy in NGC 3665. 
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Table 3. Average properties of each subsample of galaxies 


CO Type 

Ngals 

Mean 

Std. Deviation 

All 13 CO 

72 

0.111 

0.046 

All Fieldt 

62 

0.096 

0.034 

Spirals 

32 

0.098 

0.029 

Mergers 

10 

0.086 

0.020 

Field galaxies (/eo/Zioo < 0.5) 

32 

0.112 

0.036 

All ETGs 

27 

0.128 

0.081 

Interacting ETGs 

8 

0.068 

0.026 

Settled ETGs 

19 

0.154 

0.074 

Virgo ETGs 

7 

0.221 

0.064 

Eield ETGs 

20 

0.096 

0.052 

Field ETGs (/so//too < 0.5) 

8 

0.120 

0.066 


t The field sample includes all galaxies not in the Virgo cluster 

its detection mask channel by channel, to obtain an integrated spec¬ 
trum (shown in Fig. O and thereafter an line flu)0. (tabulated in 
Table[T) by summing the channels with emission (shaded blue in 
Figurej^. The un-shaded portion of the spectra in Fig.j^are 
channels where the number of pixels in the mask fall to fewer 
than the number of pixels in the beam. In these channels, the chan¬ 
nel flux is calculated by summing pixels in a contiguous region of 
the channel map the size of the beam area, and can be considered 
the baseline scatter for the stacked spectra. 

The RMS per channel was determined by calculating the stan¬ 
dard deviation of the line flux in an area equivalent to the mask, 
created by the MIRIAD task restor in “residual” mode, and mul¬ 
tiplying by VN, where N is the total size of the CO mask area 
normalized to the beam area. The RMS on the total line flux was 
then calculated by multiplying the noise per channel by the width 
of the channel At/ (shown in Table[Tll and multiplying by the square 
root of the sum of the noise squared in each individual channel (the 
shaded channels in Fig.j^: 


where i is the channel number. 

There is a further 20% absolute flux calibration uncertainty 
associated with millimeter observations, due to the time variability 
of millimeter flux calibration sources. While this 20% uncertainty 
impacts the total line flux measurements, it does not factor 

into ,^ 13 , given that all and observations were taken si¬ 
multaneously. A further 30% up-correction was made to the RMS 
to attempt to account for the oversampling of our interferometric 
maps. The 30% correction was derived by calculating the noise in 
the oversampled ( 1 " pixel) maps, and comparing that noise to the 
RMS derived from maps whose pixel area equals that of the syn¬ 
thesized beam. The RMS per channel for each galaxy is shown in 
the upper-left corner of each panel. 

Table [T] lists the galaxies that were observed in during 

the CARMA ATLAS®° survey, with the derived integrated flux 

® Line flux is defined as being the emission detected by the interferometer 
in the masking aperture, summed over the channel range determined 

by the observer 


densities and .^13 values from the stacked spectra. Of the 17 galax¬ 
ies, 15 are detected using the a priori masking technique. Figure 
presents the stacked spectra of the 15 detected galaxies. 

Comparing the line shapes of (blue) and (red), we see 
that they are generally in good agreement, supporting our claim that 
the stacking technique as well as our chosen mask threshold have 
succ essfully detected emission. 

ICrocker et al.l ( 12012h published IRAM 30m single-dish spec¬ 
tra of i^CO(l-O) and ^^ CO(2-l) for 18 of th e ATLAS®° galaxies. 
Our sample overlaps the lCrocker et al.l ( 12012h sample in the case of 
8 galaxies: IC676, NGC2764, NGC3607, NGC3665, NGC4694, 
NGC 4710, NGC 5866 and PGC 058114. We detect all overlapping 
galaxies at a signal-to-noise ratio > 3. CARMA recovers more 
emission than the IRAM 30m in all cases except IC0676, equal 
within the uncertainty (see Table [T}. In many cases, this is likely 
because the ^®CO emission subtends a larger angular size than the 
30 m beam (as was oft en the case for the ^^CO emission published 
in lAlatalo et al.ll2013h . Because ^®CO and ^^CO were measured 
simultaneously with CARMA, ,^13 does not suffer from flux cali¬ 
bration uncertainties. We therefore elect to use the ^®CO integrated 
flux densities from the current work for all overlapping galaxies. 

Radio continuum pollution leading to enhanced ^®CO fluxes 
is not likely to be an issue. NGC 3665 and NGC 4710 are the only 
two galaxies in this sample with detected continuum, at « 9 and 
4 mJy, respectively, at 110 GHz. The ^^C O data processing in cludes 
a continuum subtraction step (detailed in lAlatalo et al.ll20l^ . 

Table [^compares i^i 3 in the nuclear and extended regions of 
NGC 3665, 4429, 4435, 4710, 4753 and 5866 (see Section[3T}. In 
the first four, ,^13 is larger in both the nuclear and extended region 
than in the stacked spectrum, and vice versa for NGC 4753 and 
NGC 5866. The discrepancies between stacking and resolved ,^13 
measurements can be explained by many effects. Stacking is able 
to detect fainter ^^CO emission from a larger area than the resolved 
data. Compared to well-resolved objects, where the brightest ^^CO 
peaks were able to fill the beam, the average ,^13 from stacking 
would be lower than from the resolved maps. On the other hand, in 
cases where CARMA has marginally detected ^^CO in the maps, 
stacking would be able to increase the total ^^CO flux detected, 
resulting in a higher ,^13 for stacking than mapping. Deeper ^^CO 
observations might be able to fully explain this discrepancy, but are 
beyond the scope of this paper. 

4 DISCUSSION 

4.1 Distribution of ^®CO in the galaxies 

Figurel^shows the spatial distribution of the ^®CO emission, com¬ 
pared to the ^^CO in ATLAS^° galaxies. Only the galaxies with 
the brightest ^®CO are detectable in a spatially-resolved way. There 
are many features worthy of notice in the ^^CO distributions. 

NGC 4710 and NGC 5866 have strong, edge-on bars whose 
ki nematic features w ere detected in the ^^CO emission (see §4.2.2 
in lAlatalo et alj20T3) . In these cases, the ^®CO preferentially traces 
the bar region, particularly a nuclear disc or ring and an inner 
ring, the locations of the largest concentrations of molecular gas 
in barred systems dAthanassoula & Bureaulll99^ . NGC 4710 and 
NGC 5866 were a lso among the handfu l of galaxies detected in 
OH absorption by iMcBride et al.l ( l2015h . confirming high optical 
depths. The ^^ CO detections at the disc edges further confir m the 
conclusions of lAlatalo et al.l ( 12013h and IPavis et ^ ( l2013ll . that 
these systems are bar-(-ring systems, as the ^^CO there is likely 


^{Fi3co) = 
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Figure 4. Resolved S?,\o, in the six strong detections amongst the 

ATLAS®^ galaxies, separating the nucleus (points) and the surrounding 
gas (open circle). These numbers are likely upper limits for integrated .^13, 
as only integrated flux densities exceeding our detection threshold 

are represented. 


tracing the portions of the edge-on rings with higher column den¬ 
sities (edge brightening). Topal et al. (in preparation) conhrm this 
and study the spatially-resolved bar signatures using a multitude of 
high-density molecular tracers in both galaxies. 

The Virgo galaxies NGC 4429 and NGC 4435 are robustly de¬ 
tected in throughout the extent of the detections. In 

both cases, faithfullly traces ^^CO. NGC 4753, known to be 
at the centre of a developing group dSteiman-Cameron et al.ll992h . 
also seems to have that traces the brightest portions of '^^CO, 
though it is likely that beyond these points the emission sim¬ 
ply drops below the detection threshold (rather than being absent). 


4.2 Radial gradients in ^13 

Figure |4] show .^13 in nuclear and surrounding gas in the 6 galax¬ 
ies mapped, with the points representing the nuclear pixels and 
the open circles representing the extended emission (with resolved 
.^13 values shown in Table |2j. As discussed in 1 13.21 ratios mea¬ 
sured from the images differ slightly from those measured with the 
stacked spectra, likely due to the images missing faint emis¬ 
sion picked up in the stacked spectra. Within the selected pixels, 
these ratios are thus likely reasonable, but they are not necessarily 
representative of the entire galaxy, where the stacking method is 
better able to recover e mission below the loca l detection threshold. 

In spiral galaxies, IPaglione et alj ( 1200 ih found a .^13 gradi¬ 
ent, increasing from the nuclear region to the outlying regions on 
average, though this result had large scatter. Molecu lar discs in 
ETGs are much less extended than those in spirals dPavis et al.l 
l2013h in absolute terms, so it is possible that a .^13 variation will 
be much less pronounced in ETG discs, either due to resolution 
or because the gas is better mixed (i.e. shorter dynamical times). 
In the 6 spatially-resolved CARMA detections (shown in Fig. 0 , 
there is no decipherable trend, and in 5 of the 6 objects the nu¬ 
clear and non-nuclear .^13 values are consistent within the errors. 
NGC 3665 is the single object in which the trend is larger than the 
uncertainties, but follows the opposite trend from spirals, and could 
be due to sensitivity limitations in detecting fainter flux. 


^^CO{1-0) in ETGs 1 

There might be a slight increase in .^13 from the nuclear region 
to the off-nuclear regions in NGC 4429 and NGC 5866, but .^13 
remains flat in NGC 4710 and decreases in NGC 3665, NGC 4435 
and NGC 4753. It is however possible that we are missing a non- 
negligible amount of faint extended emission, due to the sensitiv¬ 
ity limitations of the CARMA sample. As mentioned above (i l3.2l (. 
the disagreement in .^13 between the image measurements and the 
stacked measurements is evidence that this might be the case. 

We therefore do not see any evidence for an increase in .^13 
across galactic discs in this subset of ETGs, in contrast to what has 
been seen in spirals, but we cannot rule out that such a gradient 
could exist, given the sensitivity as well as resolution limitations of 
our observations. Deeper, higher resolution follow-up observations 
should more robustly establish whether ETGs indeed lack a .^13 
gradient, whether this is an environmental effect, or whether our 
result is due to sensitivity issues. 


4.3 Variations of ^13 based on galaxy type and environment 

Figure [ 5 ] presents the distributions of .^13 as a function of galaxy 
types, in t he new CARMA de tections of Virgo and non-Virgo 
Crocker et ^ l l2012h ETGs, NGC 404 and NGC 5195 
as well as literature s p irals and merger s 


ETGs, the 

Crocker et alJ 

from ISage 

1 1990 I). as 


Figurej^also shows the .^13 distribution as a function of environ¬ 
ment in the case of ETGs only (as there is no Virgo spiral with 
^^CO data currently in the literature). Spirals and mergers tend to 
have .^i 3 < 0.1, while ETGs are distributed much more evenly as 
a function of .^ 13 . _ 

W e use the morphological classifications from I Alatalo et al.] 
( l2013h to differentiate between settled and disturbed gas morpholo¬ 
gies, though NGC 1266 was re-classified as disturbed, since the tur¬ 
bulent nature of the molecular gas is visi ble in the ^^CO and CS(2 - 
1) lines using higher resolution imaging dAlatalo et al.l201 ill2015h. 
The ^^CO in NGC 404 was classified as a ring dWiklind & Henkel 
Il990h . and is therefore settle d, while that in NG C 5195 would be 
classified as mildly disrupted dKohno et al.l2002h . and is considered 
disturbed. The .^13 value within ETGs seems to have a strong de¬ 
pendence on the disruption of the molecular gas, with disrupted gas 
distributions tending toward the lowest .^13 When ETGs are split by 
environment, the Virgo ETGs represent the majority of ETGs with 
.^i 3 > 0.1, while field ETGs average less than .^i 3 = 0 . 1 . The dis¬ 
tributions not only show that environment seems to play a role, but 
also that disrupted molecular gas distributions tend toward lower 
.^13 in ETGs. 

Eigure[^compares .^13 wi th the dust temperature, r epresented 
by the 60-to-100 )rm flux ratio. lYoung & Sanders! ( Il986h discussed 
the dependence of .^13 on the dust temperature, with lower values 
of corresponding to hig her values of .feo/.fi oo, further con¬ 
firmed by the o bservations o f Sage & Isbell (Il99 ih who increased 


the sample size. lAalto et ^ ( Il995h drastically increased the 
her of spirals for which ^®CO observations exist, confirming the 
trend as well as hypothesizing that turbulence within the molecu¬ 
lar gas might influence .^ 13 . Warmer dust temperatures are a rea¬ 
sonable tracer of turbulence in the underlying molecular gas, as an 
abundance of young stars will provide both UV p hotons to heat th e 
dust and momentum to stir up the molecular gas. dOgle et alj2010l) . 
However, the presence of Virgo galaxies at low /eo //too also hav¬ 
ing the highest indicate that environment could also be a factor 
in se tting .^ 13 . 

ICrocker et al.l ( 1201 2h found evidence that Virgo galaxies and 
those galaxies in groups have boosted values, but due to small 
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Figure 5. (Top): Distribution of all galaxies with (gray), compared with the distribution of spirals (orange), mergers (red) and ETGs (turquoise) only. Both 
spirals and mergers seem to have an average 0.1, whereas ETGs have a much flatter distribution. (Bottom): Distribution of all ETGs with .^isfgray), 

compared to the distribution of CO-disrupted ETGs (pink), field ETGs (blue), and Virgo ETGs (green) only. Disrupted ETGs tend to have the smallest 
values, while field ETGs are consistent with spirals. Virgo ETGs account for the majority of large .S’13. 


number statistics were not able to confirm tha t result. The CARMA 
detections add many more ETGs to the lCrocker et al.l l l2012h 
work, including an additional 3 Virgo galaxies. The addition of 
these ETGs not only strengthens the trend, but appears to also fur- 
ther support the suggest i on that turbulence might influence .^13 
jYoung & Sanders! 1 19861 : lAalto et aPIl995h . Magenta outlines in 
Fig. | 6 ] correspond to ETGs with disturbed molecular gas. The 
disturbed class of ETGs almost exclusively appears at values of 
.^i 3 < 0.1, the only exception being NGC4694 that also is a mem¬ 
ber of the Virgo cluster. 

In spirals and mergers, there appears to be a ceiling in .^13 
of around 0.15 for the majority of sources. Once the ETGs are in¬ 
cluded, a large proportion of them inhabit the Miz> 0.15 region, 
with Virgo ETGs exclusively lyin g above this threshold. N GC 4710 
(identified as a Virgo member in ICappellari et al.1l20l lah is found 
to have the lowest .^13 in the Virgo ETG sample, of 0.15, though 
it is of note that while NGC 4710 has a velocity similar to other 
Virgo cluster members, it is not formally within the virial radius 
of Virgo. .^13 in the CARMA sample is on average larger than 
that in all other samples discussed, due in part to the fact that the 
CARMA observations include more ETGs with predominantly set¬ 
tled gas distributions (which were not among the brightest ^^CO 
detec tions in ATLAS^° fi rst followed up in ^®CO with the IRAM 
30m: ICrocker et ai]|2012h . This inclusion of the ETGs with settled 
molecular gas seems have increased the population of objects in the 
.^i 3 > 0.15 region by over a factor of two. 

Table [Ogives the .^13 mean and standard deviation of each 
subsample of galaxies shown in Fig. |5] including isolated galaxies 
based on environment, morphology and dust temperature. Overall, 


Virgo ETGs have a larger mean .^^l 3 value than all other subsam¬ 
ples and are distinct (means inconsistent within the standard devia¬ 
tions) from all other subsamples except for cold dust field ETGs (of 
which there are only 8 ). While these two populations agree within 
the errors, the cold dust field ETG population is much more similar 
to the complete cold dust field sample than to Virgo ETGs. Thus, 
while we cannot completely rule out that the variation in .^13 that 
we observe is solely dependent on molecular gas turbulence (as 
traced by dust temperature), it does appear that ^13 in Virgo is 
boosted. 

Figure [^isolates the ETGs and spirals in the ATLAS®® par¬ 
ent sample Icappellari et alj 1201 lal) . performing S pearman rank 
statist ics as was done for the original sample from Icrocker et all 
( l2012h . The Spearman rank correlation coefficient and the prob¬ 
ability that the null hypothesis (no correlation) is true are given. 
Low probabilities for the null hypothesis (we define ‘low’ to be 
less than 0.03) suggest the data are corr elated. In most cases , the 
correlations that were originally seen in ICrocker et alj l l2012t) (in¬ 
cluding those with (1) dust temperature, (2) stellar mass and (3) 
environment) persist within the larger sample, and a new correla- 
tion is observed w ith the (4) molecular gas extent (i?moi/7?e from 
IPavis et^l2013h . which was only performed on the ETGs. The 
dependence o n S 3 (the 3 Mpc numbe r density; a proxy for envi¬ 
ronment from Icappellari et al.ll 2011 ^ has strengthened with the 
increased number statistics, but so has the dependence on galaxy 
mass. We are thus unable to determine whether the variation in 
,^13 is independently correlated with stellar mass, environment, or 
molecular extent, or if our ^®CO-detected Virgo ETGs are simply 
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Figure 6. Integrated ^^CO/'^^CO ratio for ETGs (green stars) are shown as a function of the 60-to-100/rm flux ratio (/eo/Zloo), and compared to spiral 

S ' ; s (spirals), and merge rs (c lovers). Upper limi ts ai'e shown as downwar d facing arrows. Literature spiral g alaxies and mergers are from lYoung ^_Sandersl 
. ISage & Isbelll il99lh . and lAalto et aljil995h . with ETGs (green) fromISagdfl990h . lCrocker et al.l2012r and this work. The two ETGs from IS ageli 199^ 
are NGC 5195 (the minor interaction companion to M 51 jKohno et alj2003) and NGC 404. The black dotted horizontal line shows the location of .®i3=0.15. 
Virgo ETGs are represe nted by green stars w ith a circular black background. All spirals are field spirals, as there is no Virgo spiral with ^13 literature data. 
While a few ETGs from ICrocker et all j2012h have elevated M 13 ratios, the majority have ratios similar to those of spirals and mergers from previous work. 
Moreover, these ETGs also follow the trend with dust temperature /eo/Zioo. with the highest .^13 corresponding to t he coldest dust temp eratures. ETGs 
(stars) with a magenta outline represent those that are morphologically classified as either mildly or strongly disrupted in lAlatalo et alj i2013ll . It is clear that 
ETGs with this classification are much more likely to have a low .^13, except in the case of NGC 4694, which is also a Virgo galaxy. While two of the high 
dust temperature and low .^13 ETGs from our study contain an AGN (NGC 3607 and PGC 029321), PGC058114, which also has a high dust temperature (but 
also a high .^13), does not appear to (Nyland et al., in preparation). 


more massive and compact than the ^^CO-detected field ETGs (dis¬ 
cussed in or spirals, causing the correlation. 


When we differentiate between types of molecular gas in 
ETGs, either using the CO morphological classifications from 
lAlatalo et alj ( l2013l) . molecular gas that is disturbed morphologi¬ 
cally preferentially has a much lower ,^13 than that of settled or 
aligned discs. In fact, on average morphologically-settled molecu¬ 
lar discs have ^13 values almost two times higher than those of dis¬ 
turbed molecular distributions. When we compared the ^13 to the 
kinematic misalig nment between the gas and the stars 
IPavis et aDl 2011 ). the probability of a null hypothesis was 24%, 
meaning that it appears that the disruption of the molecular gas mat¬ 
ters far more than whether the gas was externally accreted. Figure|^ 
ou tlines in magenta ET Gs that have been classified as “disrupted” 
by lAlatalo et al.l ( l2013h . and shows that in general those galaxies 
have lower ,^ 13 , (in addition to high integrated flux densi- 

ties, given they wer e part of the flux-limited dense gas survey from 
ICrocker et all^oT^l . 


This fits quite well with the work of lAalto et al.l l ll995t) . who 
see similarly low ,^13 values in mergers, compared to normal spi¬ 
rals. It is difficult to parse whether the differences in distributions 
we see are because galaxies in Virgo tend to have much more set¬ 
tled molecular discs, or if there are other effects at play. Further 


studies of in a large number of ETGs with settled molecular 
discs, in a wide variety of environments, will likely shed light on 
the relative importance of morphology and environment. No previ¬ 
ous study includes observations of Virgo spirals. 


4.4 Possible causes of the variation in .^ 13 , optical depth or 
abundance? 

In the giant molecular clouds (GMCs) of our Milky Way, typical 
. ^13 values are 0 .2 for whole GMCs, and 0.3 for GMC centres 
jPolk et alJll988h . So an ensemble dominated by cold GMC-like 
clouds should have a high ,^^l 3 . Ratios above this can be under¬ 
stood as arising from optically thick and very cold gas, perhaps like 
a GMC but with its surrounding envelope stripped away. We know 
that the ^®CO/^^CO li ne ratio correlates strongly with the dus t 
colour (or temperature: I Young & Sanderill986l : lAalto et alJll995h . 
as is clearly seen in Figure]^ Alternatively, a high ^13 value could 
be partially a result of enhanced '^^C over ^^C, although this would 
still require very optically thick gas. An illustration of the latter is 
the fact that the i^c/^^C abundance ratio d ecreases strongly with 
radius in the M ilky Wav jMilam et al.ll20 0l). but .^13 has a strong 
oppos ite trend Joigel et ^ 1990l : CMpenter et aDll990l : ISodrosld 
I 199 ll) . This is because the ^^CO optical depth in Galactic centre 
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Figure 7. ^13 against ISM properties and environment of the ^^CO-detected ET Gs from ICrocker et all i2012h and this work in addition to the spiral galaxies 
that are part of the 87 1 galaxy parent sample from which l^pnellari eta n ^lQllbl) selected the ATLAS^^ ETGs. From left to right: (1 ) /eo/ZlOO ratio of IRAS 
fluxes (a dust temperature proxy); (2) Mk (a stellar ma ss proxy); (3) S3, the 3Mpc galaxy density from lCappellari et alj i201 lal) (an environment proxy). 
(Bottom): (4) i?mol/^e, the molecular gas extent from lPavis et alJ i2013h . The Spearman rank correlation coefficient (‘r’) and the probability that the null 
hypothesis (To’, no correlation) is true are given at the top of each plot. Field ETGs are red and Virgo ETGs are blue, while spirals are green with gray error 
bars. There appears to be a correlation in each with each parameter, though the galaxies within Virgo tend to be very massive and to contain more compact 
molecular discs. Virgo galaxies also tend to have the highest ^13 even within individual mass bins. 


clouds is significantly lower due to both temperature effects and 
the presence of unbound gas. 

For the ^^CO/^^CO line ratio, optical depth effe cts may ei¬ 
ther b e caused by temperature/density changes (e.g. lAalto et ^ 
Il997l) . o r by increased lin e-widths in diffuse, non self-gravitating 
gas fe.g. lAalto et alj201(]|) . These do not require a secondary boost 
in the isotope compared to that has been suggested as tak- 
ing place during st e llar processing , speci fically in low-mass stars 
dCasoli et alJll99lh . ICrocker et alJ ( 1201 2 h considered abundance, 
fractionation and optical depth explanations for the trends seen in 
their smaller sample of ETGs, and concluded that optical depth 
(mediated by the gas dynamical state) dominates the ^13 ratio 
changes. Here, with additional sample galaxies and more informa¬ 
tion about the gas dynamics from the interferometric maps, we re¬ 
turn to this question. In addition, we consider a new option, the 
evaporation of lower-mass GMCs, or “the survival of the densest.” 


4.4.1 Fractionation 

Carbon-bearing molecules are susceptible to a process called 
fractionation. This is particularly true for CO, and therefore 
^^CO/^^CO abundance ratios may deviate from the ra¬ 

tio from nucleosynthesis. There are two competing processes. The 
deviation may be in the form of an enhancement or a reduction de¬ 
pending on which process dominates. The first process is selective 
photodissociation, which will enhance since it is the more 

abundant isotopic variant an d self-shields more effectively (e.g. 
Ivan Dishoeck & Blacklll98^ . The second process operates when 
gas kinetic temperatures are low to moderate (below 35 K) and el¬ 
evates through the isotopic charge exchange reaction: 

i2co_p i3q+—^ 1^00 +AE. 

This results in ^®CO/^^CO abundances above the 
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ratio from nucleosynthe sis by up to a factor of five (e.g. 
ISonnentrucker et ^l2007h . without also requiring an enhancement 
in the abundance. 

How can we test if our .^13 ratios are caused hy ion molecule 
exchange reactions? The dust temperatures are low, so if they are 
related to the gas kin etic temperatures fractionation seems possible. 
iRitchev et all ' 201 ih point out that because is the most abun¬ 
dant carbon-bearing molecule in the ISM, an enhancement in 
will result in a depletion of the carbon reservoir of Other 

molecules originating from the remaining carbon should therefore 
have a ratio that behaves in the opposite sense compared 

to 1 ^ 13 , i.e. we would expec t to see smaller ratios in these 

species for Virgo galaxies. IRitchev et all 1201 ih suggest CN and 
CH’*' as comparison species, and for our purposes CN should work 
well. CN has the additional advantage that the two adjacent spin 
groups can be used for an accurate measure of the optical depth. 
HCN and HCO^ should also work quite well. We note that, for ho¬ 
mogeneous clouds, the enhancement of ^®CO should occur mostly 
in the outer regions where there is sufficient ^®C^. If the clouds 
have structure the enhancement will leak further into the clouds. 

While fractionation i s a compelling possi bility to create the 
enhanced .^13 ratios seen. ICrocker et all l l 2012 h did not see a cor¬ 
relation between '^^CO/HCN (compared to '^^CO/HCN) and mor¬ 
phological type, or temperature, potential signs that fractionation 
is an important driver of the ^13 variations. Observations of dense 
gas across a larger sample of molecular gas-containing ETGs would 
further clarify whether or not we can completely rule out fraction¬ 
ation, but current data suggest it is not a dominant effect. 

4.4.2 Enhanced abundances from stellar reprocessing 

A possible origin for the distribution of .^13 ratios is that the molec¬ 
ular gas in Virgo galaxies has reside d in those galaxie s for far 
longer than that of average field ETGs. lPavis et al.l l l201lh discov¬ 
ered that the angular momentum of ionised and molecular gas in 
Virgo galaxies is always aligned with the stellar angular momen¬ 
tum, which was posited to mean that the molecular gas within Virgo 
ETG s is not externally accret ed, and either is due to stellar mass 
loss jFaber & Gallagheilll976h or remains from aligned molecular 
gas when the galaxy first entered Virgo, thus supporting this claim. 

The observation that some Virgo spirals may have signifi¬ 
cantly hi gher O/H abundance s within their ionised gas than field 
galaxies dShields et alj Il99lh also support this claim, althoug h 
more recent metallicity stu dies contradict it dHughes et alj|201^ . 
ICasoli et^ dl99lL Il992al l9l discuss low .^13 ratios as likely be¬ 
ing an abundance effect, with more primordial (and largely unpro¬ 
cessed) molecular gas containing fewer atoms, since this iso¬ 
tope is formed via low-mass stellar nucleosynthesis. The interstel¬ 
lar medium (ISM) in Virgo galaxies would have existed sufficiently 
long for this enrichment to occur. NGC4710, which is on the out¬ 
skirts of Virgo, has .^13=0.15, at the junction between most spirals, 
interacting galaxies, and Virgo galaxies. If NGC4710 is on its first 
approach (suggested by its presence outside of the Virgo virial ra¬ 
dius), and therefore just losing the ability to accrete new primordial 
gas, then this enrichment has just started to boost the abun¬ 
dance and is beginning to impact .^ 13 . 

Objects with morphological signs of recent accretion events 
are the most likely to also have low ^13 values, as is expected if 
the gas in the objects that recently accreted is more primordial and 
less enhanced. Though this is the case, kinematic misalignments 
with stars, consistent with external accretion, did not show such a 
distinctive difference, especially compared to the disruption seen 


in the gas (see il4.3b . If we assume that on average, the accreted 
molecular gas has a lower abundance, then there should be a 
trend toward lower .^13 from externally acquired gas, not solely 
morphologically disrupted gas. 

Studying isotopic line ratio in other gas species (HCN vs. 
H^^CN, for instance) might be able to confirm whether enrichment 
plays a critical role in the .^13 ratio. Alternatively, studying .^13 at 
higher excitation, and getting a different value of .^13 in the settled 
discs, could rule out enrichment (though not necessarily confirm 
it). 

4.4.3 Midplane pressure 

Rather than being an abundance effect, the elevated .^13 in Virgo 
ETGs could also be the result of midplane pressure, that can 
modify the optical depth of the molecular gas. First posited by 
ICrocker et alj l l 2012 h to describe the correlations seen in the orig¬ 
inal sample of ATLAS®® ETGs, it is possible that the hot intra¬ 
cluster medium (ICM) within Virgo produces a larger than normal 
midplane pressure in galaxies. This excess pressure on the molec¬ 
ular disc would compactly it, creating a disc with higher optical 
depth (and thus higher .^ 13 ) than is seen in field galaxies. The fact 
that the resolved maps of NGC 4429 and NGC 4435 do not appear 
to show any significant spatial variation in .^13 seems to support 
the possibility of a high pressure throughout. 

However, high optical depths within the discs of these galaxies 
would likely imply a higher frac tion of dense gas, and thus i mpact 
their s tar formation efficiencies jGao & Solomoiiir2004tl . but lPavij 
( I 2 OI 4 I) have shown that the star formation efficiencies of Virgo 
ETGs are similar to those of all ETGs. Since it also appears that 
the molecular gas distribution plays a role in this ratio for ETGs, 
the high .^13 in non-Virgo settled ETGs seems to indicate that mid¬ 
plane pressure is not the dominant effect in these systems. 

4.4.4 Preferential evaporation of lower density GMCs 

It is possible that the observed enhancement of .^^’13 is nothing more 
than “survival of the densest” (GMCs). As Virgo galaxies move 
through the ICM, and the low density clumps are stripped (leaving 
only the densest clumps, likely with high .^ 13 ), one would expect 
our .^13 measurements to prob e only the hi g hest e nd of the usual 
density distribution of GMCs. IPineda et"^ (l2008l) show that the 
highest .^13 values in local clumps (.^13 « 0.4), with derived den¬ 
sities of 10 ® — lO'* cm~®, are indeed comparable to the average 
.^13 of the Virgo gala xies. The mild correla tion of the HCN/^^CO 
ratio with E 3 seen by ICrocker et alj ( 1201 2 h also supports this hy- 

i iothesis. In fact, r esults from the Herschel Virgo Cluster Survey 
Davies et al .l2012h suggest that dust and gas densities are higher in 
Virgo cluster galaxies than in field galaxies, likely requiring denser 
constituent molecular clouds. Higher resolution studies of GMCs 
in Virgo galaxies with ALMA would further elucidate whether the 
properties of the dust transfer directly to the molecular gas proper¬ 
ties in these systems. 

4.5 Investigating the origins of the .^13 variation 

To fully understand the dependence of .^13 on environment, as well 
as what mechanism drives that dependence, further ^®CO obser¬ 
vations of Virgo galaxies, including Virgo spirals, will likely be 
required, hopefully shedding light on what is causing the large 
spread in .^^’13 values observed (including the Virgo-field galaxy 
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dichotomy). If abundance is the most important factor, then 
the gas phase metallicity (a proxy for “primordality”) should cor¬ 
relate with ,^ 13 . If midplane pressure is the dominant factor for 
the variation in then the ratio in varying environments 
(including galaxy groups and possibly the Coma Cluster) should 
correlate. Finally, if entry into the hot X-ray gas halo of the Virgo 
cluster has stripped all but the densest clumps, then a measure of 
how long each galaxy has been within the cluster should correlate 
with ,^ 13 . A high reso lution census of the GMCs in Virgo (similar 
to the investigation of lUtomo et al.ll2oT5h might also provide evi¬ 
dence for whether “survival of the densest” is able to explain the 
variations we see in ,^ 13 . 


5 CONCLUSIONS 

The CARMA ATLAS®° ^^CO survey of I Alatalo et Zl (1201 3l) in¬ 
cluded 17 galaxies simultaneously observed in ^^CO(l-O), increas¬ 
ing the number of ATLAS^° ETGs observed in ^^CO by 50%, and 
including far fainter ^^CO detections (and mo re objects w i th set ¬ 
tled molecular gas) than the original work of iKrips et al.l l l2010ti 
an d ICrocker et al.l l l 2012 h . 

We presented here the ^®CO(1-0) resolved maps of 6 ETGs 
from the CARMA ATLAS^° survey, those detected with sufficient 
significance to derive spatial-resolved ^®CO i nformation. An a pri - 
ori masking analysis based on the method of lSchruba et al.l l l201lh 
allowed us to detect ^®CO in 15/17 of the observed sources, this in 
a spatially-integrated manner. 

There does not appear to be evidence for spatial variations in 
,^13 in the 6 ETGs we were able to map, though it is likely we 
would not hav e been able to detect a gradient at the level found in 
IPaglione et all ( 1200 ih . given the signal-to-noise and resolution of 
our data. 

Combini n g the se new observations with those from 
ICrocker et alj d 2012 h and other data from the literature, we 
find that Virgo ETGs have a significantly higher .^13 than both 
field ETGs, and spirals and interacting galaxies, although ETGs 
with settled molecular discs also have preferentially higher ,^13 
than both spirals and ETGs with morphological disruptions. It 
is also possible that the variations of .^13 are driven by mass or 
molecular extent, as Virgo ETGs are, on average, more massive 
and contain more compact molecular discs than field ETGs. 

We hypothesize that the higher .^13 ratio in Virgo could be 
caused by extra low-mass stellar enrichment taking place in Virgo 
cluster galaxies (boosting the ^®C elemental abundance in the ab¬ 
sence of external gas accretion), by increased midplane pressure 
due to the ICM, or by the survival of only the densest clumps of 
molecular gas (as the galaxies enter the Virgo ICM). Studies of 
^^CO in more Virgo galaxies, particularly spirals, as well as studies 
of other ^®C and ^^C isotopologues, should be able to distinguish 
between these hypotheses. 
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